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ABSTRACT 

We present a search for radio transients in tlie field of tlie briglit radio source 
3C 286 using arcliival observations from tlie Very Large Array. These observations 
span 23 years and include 1852 epochs at 1.4 GHz in the C and D configurations. 
We find no transients in the field. The sensitivity of the observations is limited 
by dynamic range effects in the images. At large flux densities (> 0.2 Jy), single 
epoch observations provide a strong limit on the transient surface density. At 
flux densities near the dynamic range threshold, we use the requirement that 
transient sources must appear in consecutive epochs to be confirmed as real. 
This sets the sensitivity at low flux densities to transient durations of r ~ 1 day, 
while r > 1 minute for high flux densities. At 70 mJy, we find a 1-a limit on the 
surface density S < 3 x 10"^ deg~^. At 3 Jy, we find a 1-a limit S < 9 x 10""^ 
deg~^. A future systematic search of the VLA archives can provide one to two 
orders of magnitude more sensitivity to radio transients. 

Subject headings: radio continuum: general — radio continuum: stars — radio 
continuum: galaxies — surveys 

1. Introduction 

Radio transient (RT) sources probe the high energy population of the Universe. Known 
hosts to transient radio emission include neutron stars, black holes, supernovae, gamma-ray 
bursts, and highly magnetized stars and planets. Most of what is known about RTs has 
been learned from follow-up of event s discovered at high-en ergy or optical wavelengths or 



through serendipitous discovery (e.g.. iBrunthaler et al.ll2009l ). Systematic searches for RTs 



have been conducted but vast parameter space remains unexplored. 
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Blind searches for RTs are an important scientific goal for maior new radio tel escope 
facilit ies such as the Al len Telesco pe Array (ICroft et al.l l2010bl : iBower et al.l |2010[ ). LO- 



FAR (Hesse: 



( iLazio et al. 



s et al.ll2009h . ASKAP flMacquart et al 



201o[ ). and the Long Wavelength Array 



2010l ). But t here is signifi c ant o pportunity to identify RTs through analy- 



sis of archival radio data. iBower et al.l (j2007| ) used nearly 1000 observations of a blank 
field observed by the Very Large Array (VLA) over 20 years t o identify a set o f RTs that 
have no apparent counterpart at radio or optical wavelengths. lOfek et al.l (120101) have sug- 
gested that these may be due to neutron stars. Levinson et al. ( 20021 ) arid Gal- Yam et al. 



( 120061 ) conducted a search for RTs at 1.4 GHz through a comparison of the VLA NVSS and 
FIRST survey catalog s that identified a radio supernova in the nearby galaxy NGC 4216. 
Bannister et al.l (120101 ) recently completed a search at 843 MHz of the Molongolo Observa- 
tory Synthesis Telescope (MOST) data archives that uncovered 15 RTs and a larger number 
of variable source s. Some of the sour ces discovered in the Molongolo search appear similar to 
the RTs fourid bvlBower et al.l (120071 ) in that they have no faint radio or optical counterpart. 
Becker et al.l (120101 ) found a population of faint, variable radio sources in the galactic plane, 
the majority of which are without multi-wavelength counterparts. 

We describe here a search for radio transients in the field of the quasar 3C 286 with 
analysis of 1.4 GHz archival data from the VLA. One of the standard fiux calibrators for the 
VLA, 3C 286 has beei i observed th o usand s of times over the life time of the array. This search 
builds on the work of iBower et al.l ( 120071 ): however, observations in the vicinity of the bright 
(15 Jy) source 3C 286 place limits on the sensitivity that can be achieved. The dynamic 
range of these observations is limited by calibration errors and other systematic effects rather 
than by statistical noise. These calibration and systematic errors can produce an apparent 
source in the image that is well above the theoretical detection threshold. Accordingly, we 
require stronger evidence (such as appearance in consecutive epochs and comparisons with 
other cataloged sources in the field) to demonstrate the existence of an RT. We present the 
data and its analysis in §2, a simulation demonstrating the ability to identify sources with 
these methods in §3, our source detections and transient identification efforts in §4, limits in 
transient surface density in §5, and a summary in §6. 



2. Data and Analysis 

All of the data used were archived VLA observations. We limited our search to 1.4-GHz 
continuum-mode observations in C and D configurations. Continuum mode has a bandwidth 
of 50 MHz in two separate intermediate frequency bands and in both circular polarizations. 
In the longer baseline VLA configurations (A and B), bandwidth smearing reduces the area 
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that can be imaged without significant losses. We selected all observations between April 
1984 and May 2007. After culling epochs with corrupt or too few data, there were 1852 
observations to inspect. The median duration of the epochs is 2 minutes; 0.2% of epochs 
have a duration shorter than 20 seconds and 0.7% of epochs have a duration longer than 20 
minutes. The total integration time included is 167 hours. 

We used an automated fiagging, calibrating and imaging pipehne ^y ithin the AIPS pack- 



age originally developed for another archival field ( iBower et al.l 120071 ). We self-calibrated 
each epoch assuming a point source at the field center. As demonstrated in §3, this method 
permits recovery of other sources over a wide range of fiux densities. Each epoch was then 
imaged over a square area slightly larger than two times the primary beam diameter. The 
presence of a bright calibrator (3C 286) at the center of the field permits us to perform 
amplitude and phase self-calibration on the data. 

The median rms fiux density of individual epoch images is 2 mJy, corresponding to a 
dynamic range of ~ 10^ (Fig.HD. We also plot the histo gram of separation between sequential 
epochs. The characteristic separation between epochs is ~ 1 day. Typical resolution in the 
C and D configurations is ~ 15 arcsec and 45 arcsec, respectively. 

We created a deep image (Fig. [2]) by combining 1517 of the epochs with the MIRIAD 
software package. The rms of the deep image excluding bright sources is 0.8 mJy and has a 
synthesized beam of 22 x 16 arcsec^. The image is dynamic range limited and probably could 
be improved with a significant self-calibration effort. However, the image is very effective at 
identifying sources in the field that may show up in individual epoch images. We identify 
10 sources in the field including 3C 286 (Table [1]). As discussed below, we tabulate the 
detections of each of these sources in all of the epochs. We report the mean fiux density 
over all the epochs, the modulation (i.e., root mean square variation) of that fiux density, 
and the number of epochs in which sources are detected (Ndet)- For sources with N^et = 0, 
we report the fiux density from the deep image as the mean fiux. The overdensity of strong 
sources in the vicinity of 3C 286 gives a useful benchmark against systematic errors. 



3. Simulation 

A concern in analysis of data of this kind is that self-calibration on the bright point 
source in the field will distort or disappear a transient source that is not included in the 
self-calibration model. The persistence of the other sources in the field demonstrates that 
self-calibration is robust against this problem for sources with fiuxes between 30 and 200 
mJy. We performed a simple simulation to explore the effect on sources with flux densities 
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from 10 mJy to 25 Jy. The simulation creates fake data using the MIRIAD task UVGEN for 
the VLA D configuration of a 2-minute 1.4 GHz snapshot observation of 3C 286 at transit. 
The model sky is composed of the deep sources in Table [T] with an additional source at a 
separation of 700 arcsec from the field center with flux Sfake- Self-calibration is performed 
with a simple model of a point source at the phase center. The data are then imaged, 
cleaned, and restored. The flux density of the source at the phase center is consistent with 
the 3C 286 flux density until Sjake ~ S^c2?,6- The fake source is recovered consistently at 
the input flux density for Sjake < 10 Jy. This simulation confirms that we can recover RTs 
in the field over a wide range of flux densities. 



4. RT Identification 

The AIPS task SAD (Search And Destroy) was used to identify all sources in the fields 
brighter than 5a over the full area imaged. We identified a total of 30067 sources in all 
epochs. We perform several exclusions to identify any transient sources. We reject any 
sources that are in the outer 2.5% of the image which may be affected by edge problems in 
the image. We reject any sources outside of the two times the primary beam radius. We 
reject any sources with a fitted size larger than 120 arcsec since we are only interested in 
point-like sources. We reject any sources within 0.08 deg of 3C 286. The latter step removes 
a large number of sidelobes associated with the bright point source. We also remove any 
sources that are matched (within a radius of 30") to the known steady sources in Table [H 
we tabulate the number of 7a detections for each of these sources, Ndet- Three sources are 
detected in nearly all of the epochs; the faintest has a flux density of 73 mJy. Variations in 
the flux density indicate systematic uncertainty of ~ 5%. None of the sources show strong 
intrinsic variability. 

After these cuts, we are left with 608 sources in the catalog. 186 of these sources are 
found to be repeating in different epochs and 3 are found to be repeating in consecutive 
epochs. Examination of the images reveals that many of these sources are associated with 
sidelobes of the synthesized beam. If we apply a higher threshold of 7a, then we reduce the 
number of single epoch candidates to 78. Only 7 candidates are seen to repeat in any epoch 
and none are repeating in consecutive epochs. The repeating candidates all cluster near the 
brightest of the steady sources (J133148+303148) and are likely to be sidelobes. 

Assuming purely Gaussian noise, the expected number of false positive sources given a 
statistical threshold a for all epochs is Njp = 1/2 erfc{a/\/2)Ntry, where erfc is the comple- 
mentary error function. Nfry = Nepoch^f/^b is the product of the number of epochs with the 
ratio of the search area to the synthesized beam area. For this experiment, Ntry ~ 6 x 10^. 
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For cr = 5, we expect Nfp ^ 20. For cr = 7, Nfp ^ 8 x 10"^. Thus, the remaining Tex 
candidates are either real RTs or they are systematic errors. Estabhshing the reahty of any 
single epoch candidate is very difficult given the variable systematic errors present in the 
data. Thus, at the low sensitivity end, we primarily make use of only the consecutive event 
counts to estimate transient sensitivity. This has a characteristic time scale of the epoch 
separation ~ 1 day. At higher sensitivity, the characteristic RT time scale sampled is ~ 1 
minute. 

We note one unusual candidate source that showed up in a single epoch (13 July 1996). 
The source appeared as a bright (422 mJy) source with point-hke structure. The steady 
sources were present in the image at the right flux densities. This candidate had the appear- 
ance of a convincing RT. However, examination of the visibility data during this 2-minute 
observation indicated a glitch. It appears that the online flagging system was effective in 
removing most of the bad data but some bad data remained. In those bad data, the phase 
center of the image was shifted. Imaging all of the data together led to point sources at the 
position of 3C 286 and at the shifted phase center. The result was an apparently convincing 
RT candidate until we removed two 10- second integration s. Note that this glitch differs from 



the end of record problem identified by lOfek et al.l (120101 ). After exclusion of this event and 



four other single epoch candidates in which imaging revealed poorly flagged and/or poorly 
calibrated data, there were no single epoch candidates brighter than 250 mJy. 



5. Transient Rate Estimates 



We calculate the transient surface density, S, which is equivalent to the two-epoch 



(Bower et al. 


2007; 


Croft et al. 


2010b: 


Bower et al. 



20101 ). The key step is to estimate the area per epoch as a function of flux density threshold. 
We can estimate this based on the image rms and the imaged area with a sensitivity above 
the threshold; this is the statistical limit given in Fig. [31 We plot here the inverse of the area, 
which is proportional to the surface density. This statistical method, however, ignores the 
effects of systematic errors caused by 3C 286 in the image. We use an alternative method of 
estimating the area available by counting the number of epochs in which the bright steady 
sources in the image are detected (Tabled]). Under the assumption that this threshold is much 
higher than the statistical noise threshold (which holds in most cases), then the total area is 
the number of detections times the field of view imaged; this is labeled as systematic limit 
for 1 epoch. We can also use the requirement of consecutive detections of the steady sources 
to estimate the area for repeating sources; this is the systematic limit for 2 epoch detections. 
The inverse area rises steeply at flux densities below 70 mJy. For flux densities greater than 
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70 mJy, we find good agreement between the statistical and systematic approaches and httle 
difference between the one and two epoch areas. The discrepancy between the statistical 
and systematic limits at low fiux densities is the penalty that we pay for using data with a 
dynamic range limit. 

We estimate the la upper limit on S and plot it against other measurments and limits 
(Figure H]). We use the systematic two-epoch limit for fiux densities < 200 mJy and then 
interpolate to the statistical limit for very large fiux densities. The surface density above 
200 mJy applies to transients with timescales > 1 minute. We find E < 3 x 10~^ deg~^ at 
70 mJy (for r ~ 1 day) and S ^ 9 x 10~^ deg~^ at 3 Jy (for r ~ 1 min). 

These limits are an improvernent oy e r limits on very bright RTs relative to the ATATS-I 
and MOST surveys dCroft et all boiObI : iBannister et aP boioh and are comparable to the 
limits from the NVSS-FIRST comparison ( Gal- Yam et al.l 2006 ). The ATATS-II results 
are a factor of ~ 3 more sensitive than the 3C 286 results at a fiux density of 300 mJy 



( ICroft et al.ll2010al ): however, the 3C 286 results are a factor of 4 times more sensitive in fiux 
densi ty. Our limits fro r n 3C 286 observations are more sensitive than the published limits 
from iMatsumura et al.l (120091 M09), however, there is considerable uncertainty about what 



is the best limit based on the M09 data. Previous limits from the same group (iKida et al. 



20081 ) also disagree with M09. We plot the stated M09 surface density from two years of 
surveying of 8.7 x 10~^ arcmin"^ = 3 x 10~^ deg~^. Our independent calculation based 
on data presented in the paper of 3 transients detected in 50 days of observing with an 
instantaneous field of view of 500 deg^ suggests a surface density of 10"'^ deg~^; using 9 
transients from 2 years of observing we estimate S ~ 2 x 10~^ deg~^. However, we cannot 
be certain that we know the duty cycle of this observing and so take the stated M09 limit of 
3 X 10^^ deg^^ as the benchmark value. For this value, we would expect to find 3 1-Jy RTs, 
which we do not find. 

The surveys compared in this plot cover an order of magnitude in frequency (0.84 to 
8.4 GHz) and several orders of magnitude in time scale that are probed (1 minute to 1 
year). Thus, each survey takes a different slice of parameter space that may probe very 
different physics and source populations. Broadly, these observations are sensitive to syn- 
chrotron phenomena, such as the explosive ejecta of radio supernova and gamma-ray burst 
afterglows, as well as X-ray binaries and active galactic nuclei. Models for most of these 
source populatio ns indicate surfac e densities that are to the left and below the dashed line 
in the plot (e.g., iRossi et al.l 120081 ). 
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6. Summciry 

We have presented an analysis of 1852 epochs of VLA observations spanning 23 years 
of the bright cahbrator 3C 286 at 1.4 GHz. This data set provides an important search for 
radio transients brighter than 70 mJy. We do not find any transients, in contradiction with 
optimistic estimates of transient surface density from M09 but consistent with hmits from 
other surveys. Differences between surveys may be a function of observing frequency, regions 
of sky covered, and systematic problems in recovering transients. 

The results are limited significantly by dynamic range of the imaging. If systematic 
errors had not contributed to the imaging, our surface density limit would apply at flux 
densities that are an order of magnitude lower. Nevertheless, the results demonstrates that 
searches around bright calibrators can provide unique information. Future searches of the 
VLA archives can improve on these results through the use of the larger number of observa- 
tions of other standard calibrators such as 3C 48 and of observations at other frequencies. 
The examination of fainter calibrators is an important way to get closer to the statistical 
noise limits under the assumption of a fixed dynamic range limit. Finally, more complete 
models used in self-calibration may permit higher dynamic range imaging. If 1% of the data 
from the VLA archives consist of calibrators suitable for transient searching, we will have 
2000 hours of usable data from the past 25 years, corresponding to an order of magnitude 
increase in sensitivity to radio transients. 

The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. 
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Fig. 1. — Histograms of survey properties. Left: Histogram of epoch RMS noise levels. 
Right: Histogram of number of days between sequential epochs. Multiple observations on 
the same day are considered to have a zero day separation while all separations greater than 
50 days are included in the rightmost bin for the sake of clarity. 
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Fig. 2. — Deep image of the 3C 286 field. 3C 286 is at the center of the field. The gray scale 
goes from 2 to 4 mJy. 
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Fig. 3. — Inverse area for statistical and systematic limits as a function of flux density 
threshold. See text for details of this calculation. 
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Fig. 4. — Transient surface density from this and other surveys as a function of flux density. 
Result from this survey is labeled 3C 286. Curves and lines indicate det ected values and 
upper limits from a deep VLA search (B07,B07. 1,B07.2: iBower et alJl2007l ). the comparison 



of the 1.4 GHz NVSS and FIRST surveys (G06: lGal-Yam_etJl20Q6[), from additional VLA 



searches (COS and F03: ICarim et al.lboosi iFrail et al.,.2003il . from the first and second ATATS 



papers fA TATS-I and ATAT S-II: ICroft et al 



2010b 



af), frorn the fir st da ta release of PiGSS 



(PiGSS-I; iBower et al.l l2010h. from the iMatsumura et al.l (M09; 120091 ) survey, and from 
the MOST search (IBannister et al.ll2010l ). The dashed line is proportional to S^^'^ and is 
normalized to B07 estimates. Lines with arrows indicate la upper limits; otherwise the 
results are indicative of detected transients (B07, MOST, and M09). 



